Dark Matter Results from 100 Live Days of XENON100 Data 
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We present results from the direct search for dark matter with the XENON100 detector, installed 
underground at the Laboratori Nazionali del Gran Sasso of INFN, Italy. XENON100 is a two-phase 
time projection chamber with a 62 kg liquid xenon target. Interaction vertex reconstruction in three 
dimensions with millimeter precision allows the selection of only the innermost 48 kg as ultra-low 
background fiducial target. In 100.9 live days of data, acquired between January and June 2010, no 
evidence for dark matter is found. Three candidate events were observed in the signal region with 
an expected background of (1.8 ±0.6) events. This leads to the most stringent limit on dark matter 
interactions today, excluding spin-independent elastic WIMP-nucleon scattering cross-sections above 
7.0 x 10~ 45 cm 2 for a WIMP mass of 50GeV/c 2 at 90% confidence level. 
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Weakly Interacting Massive Particles (WIMPs) are a 
well-motivated class of particles [1] to constitute a major 
fraction of the dark matter in the Universe [2]. These 
particles can be searched for in underground-based de- 
tectors [3] through their coherent scattering off target 
nuclei. A quasi-exponentially falling energy spectrum 
of nuclear recoils is expected, extending up to a few 
tens of keV at most for WIMPs scattering off a xenon 
target [4]. A wide range of WIMP-nucleon cross sec- 
tions predicted by theoretical models remains untested by 
current-generation searches, with the most stringent lim- 
its on the elastic spin-independent WIMP-nucleon cross- 
section coming from CDMS-II [5], EDELWEISS-II [6] 
and XENON100 [7]. 

XENON100 is the current phase of the XENON dark 
matter program, which aims to improve the sensitivity 
to dark matter interactions in liquid xenon (LXe) with 
two-phase (liquid/gas) time-projection chambers (TPCs) 
of large mass and low background. A key feature of 
XENON100 is its ability to localize events with mil- 
limeter resolution in all spatial dimensions, enabling the 



selection of a fiducial volume in which the radioactive 
background is minimized. The simultaneous detection of 
charge and light signals provides discrimination between 
the expected WIMP-induced nuclear recoil (NR) signal 
and interactions from the electromagnetic background in 
the form of electronic recoils (ERs). First results from 
XENON100 [7] have shown a sensitivity competitive with 
that from the full 612 kgxdays exposure of CDMS-II [5] 
after only 11.2 live days of data taking. Here, the result 
of a WIMP search using 100.9 live days of XENON100 
data is reported. 

XENON100 is filled with 161 kg of ultra-pure LXe. 
Of these, 99 kg are used as active scintillator veto, 
surrounding the optically separated 62 kg target in 
47T. Thanks to careful material selection [8] and detec- 
tor design, XENON100 achieves an experimentally ver- 
ified ER background in the relevant low-energy region 
of < 5 x 10~ 3 eventsx (keV ee xkgxday) _1 (keV ec = keV 
electron-equivalent [9]) before signal discrimination [10]. 
An interaction in the cylindrical LXe target of ^30 cm 
height and ~30 cm diameter generates prompt scintilla- 
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tion light (SI) and ionization electrons, the latter being 
detected through the process of proportional scintilla- 
tion (S2) in the gaseous xenon above the liquid. Both 
SI and S2 signals are registered by photomultiplier tubes 
(PMTs), at the bottom of the LXe target for optimal 
light collection, and placed above in the gas phase. The 
interaction vertex is reconstructed in 3 dimensions, with 
the (x, y)-position determined from the hit pattern of the 
localized S2 signal on the top PMT array, and the z- 
coordinate deduced from the drift time between the SI 
and S2 signals. This allows the fiducialization of the tar- 
get volume to exploit the excellent self-shielding capabil- 
ities of LXe. Due to their different ionization densities, 
ERs (7, f3 background) and NRs (WIMP signal or neu- 
tron background) have a different S2/S1 ratio, which is 
used as discrimination parameter. 

The 242 PMTs used in XENON100 are T'-square 
Hamamatsu R8520-AL PMTs with a quantum efficiency 
of ^30% at the Xe light wavelength of 178 nm, and low 
intrinsic radioactivity [8]. The measured average energy 
threshold of the LXe veto is ~ 100keV ee - 

The TPC is installed inside a vacuum insulated stain- 
less steel cryostat which is surrounded by a passive shield 
made of high purity copper, polyethylene, lead and water 
in order to suppress external backgrounds. A constant 
flow of high-purity nitrogen boil-off gas keeps the 222 Rn 
level inside the shield < lBq/m 3 . A 200 W pulse tube 
refrigerator, installed outside the shield structure, keeps 
the detector at its operating temperature of — 91°C, with 
excellent stability over time (fluctuations <0.05%). To 
bring calibration sources ( 60 Co, 137 Cs, 241 AmBe) close 
to the target, a copper tube penetrates the shield and 
winds around the cryostat. XENON100 is installed un- 
derground at the Italian Laboratori Nazionali del Gran 
Sasso (LNGS) below an average 3600 m water equivalent 
rock overburden, which reduces the muon flux by a fac- 
tor - 10 6 . 

At low energies, the event trigger is provided by the S2 
signal. The summed signal of 84 central PMTs is shaped 
and fed into a low-threshold discriminator. The trigger 
efficiency has been measured to be > 99% at 300 photo- 
electrons (PE) in S2. 

Three algorithms are used to reconstruct the (x,y) co- 
ordinates of the events. They yield consistent results out 
to a radius of 14.2 cm, with the active TPC radius be- 
ing 15.3 cm. The (x,y) resolution was measured with 
a collimated source and is <3 mm (lcr). The algorithm 
based on a Neural Network gives the most homogeneous 
response and thus is used for event positioning, while 
the information from the other algorithms is used for 
consistency checks. The drift time measurement gives a 
z-position resolution of 0.3 mm (lcr) and allows to dis- 
tinguish of two interaction vertices if separated by more 
than 3 mm in z. The positions are corrected for non- 
uniformities of the drift field, as inferred from a finite- 
element simulation and validated by data. 



XENON100 uses continuous xenon purification 
through a hot getter. The mean electron lifetime r e is 
indicative of the amount of charge lost to impurities [11]. 
It increased from 230 /is to 380 /is for the data reported 
here, as measured weekly with 137 Cs calibrations. A 
linear fit to the r e time evolution yields the ^-correction 
for the S2 signals with negligible systematic uncer- 
tainty (< 2.5%). {x,y) variations of the S2 signal are 
corrected using a map obtained with the 662 keV ee line 
from 137 Cs. 

The spatial dependence of the SI signal due to the 
non-uniform light collection is corrected for using a map 
obtained with the 40 keV C0 line from neutrons scatter- 
ing inelastically on 129 Xe. It agrees within 3% with 
maps inferred from data using the 662 keV ee line and the 
164keV ee line, from neutron-activated 131m Xe. The light 
yield L y (122keV ce ) = (2.20 ± 0.09) PE/keV cc at the ap- 
plied drift field of 530V/cm in the LXe is determined 
by a fit to the light yields measured with all available 
calibration lines [7]. 
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FIG. 1: All direct measurements of £ e & [12, 13] described by 
a Gaussian distribution to obtain the mean (solid line) and 
the uncertainty band (la and 2a). Below 3keV nr the trend 
is logarithmically extrapolated to C c s = at lkeVnr- 

The NR energy E BI is inferred from the SI signal us- 
ing E ni = (Sl/L y )(l/C e s)(S ee /S ni ). The scintillation ef- 
ficiency £ e ff of NRs relative to that of 122 keV co 7-rays 
at zero field is taken from the parametrization shown in 
Fig. 1, which is strongly supported by the most recent 
measurements by the Columbia group [12] but includes 
all direct measurements of this quantity [13]. £ e g- is log- 
arithmically extrapolated below the lowest measured en- 
ergy of 3 keV nr , motivated by the trend in the data points 
as well as studies which simultaneously take into account 
light and charge signal [14]. The electric field scintillation 
quenching factors for ERs S cc = 0.58 and NRs S nl = 0.95 
are taken from [15]. 

The dark matter data presented here was acquired be- 
tween January 13 and June 8, 2010. About 2% of the 
exposure was rejected due to variations in detector op- 
eration parameters. In addition, 18 live days of data 
taken in April were rejected due to an increased elec- 
tronic noise level. With an average data acquisition live 
time of ~ 90%, and regular calibration runs during the 
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data-taking period, this leads to a dataset of 100.9 live 
days. This data had been blinded below the 90% ER 
quantile in log 10 (S2/Sl)-space for SK160 PE. 

It was decided a-priori to derive the dark matter re- 
sult based on a Profile Likelihood analysis as introduced 
in [17] but taking into account all relevant backgrounds 
for this dataset. This analysis does not use a cut based 
on the S2/S1 discrimination. Both, the signal and the 
background hypothesis would be tested regardless of the 
observed data. In parallel, an analysis based on the op- 
timum interval method [18] was also performed. The re- 
stricted S2/S1 space used in the latter analysis, together 
with the energy interval, defines a benchmark WIMP 
search region that allows to directly compare the ob- 
served signal with the expected background. 

From a comparison of the measured background rate 
with Monte Carlo simulations of the XENON100 elec- 
tromagnetic background [10], a nat Kr concentration of 
(700 ± 100) ppt is inferred for the data reported here, 
higher than reported earlier [7]. The additional Kr was 
introduced by an air leak during maintenance work on 
the gas recirculation pump, prior to the start of the 
data-taking period. This results in an expected ER back- 
ground of < 22 x 10~ 3 evtsx(keV ee xkgxday) _1 before 
S2/S1 discrimination, more than an order of magnitude 
lower than the one of other dark matter search exper- 
iments. After the science run presented here, the Kr 
concentration in the Xe has been reduced by cryogenic 
distillation to the level reported in [7] , as confirmed with 
a /3-7-coincidence method [16]. XENON100 is taking new 
data with this reduced background and improved perfor- 
mance. 

The imposed requirements to the quality and topol- 
ogy of events are designed to retain the highest possible 
acceptance of the expected WIMP-induced single-scatter 
NRs. The majority of cuts were designed and fixed before 
unblinding the signal region, based on expected signal 
characteristics, on NR data from a calibration of 2.9 live 
days with a 241 AmBe source, and on low-energy ERs from 
Compton-scattered gammas from a 60 Co source, which 
were recorded over a total of 5.8 live days spread in time. 
To reject interactions with a very high energy deposition, 
the identified SI and S2 signals are required to amount 
for more than half of the total digitized signal. To satisfy 
the primary requirement of the WIMP signature to be a 
single scatter interaction, one S2 signal above 300 PE 
is required, corresponding to about 15 ionization elec- 
trons. All other S2 signals have to be small enough to 
be consistent with PMT afterpulsing or delayed ioniza- 
tion signals from the same single-scatter interaction. The 
corresponding SI signal must be above 4 PE and must 
satisfy a two-fold PMT coincidence in a ±20 ns window, 
without having a coincident signal in the LXe veto. Any 
other Sl-like signal must be consistent with electronic 
noise or unrelated to the S2, based on its S2/S1 ratio. In 
addition, both SI and S2 PMT hit patterns as well as the 



width of the S2 pulse are required to be consistent with 
a single interaction vertex at the reconstructed position. 
The cumulative cut acceptance, used by both analyses, 
is shown in Fig. 2 and has an error of ~ 3%. It is esti- 
mated based on Monte Carlo simulations, 241 AmBe and 
60 Co calibration data, as well as ERs recorded outside 
the WIMP search region during the dark matter search. 
It includes a WIMP mass m x dependent S2 acceptance 
which is derived from the expected recoil spectrum and 
the measured S2 vs. SI distribution. 
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FIG. 2: Acceptance of all data quality cuts used for the analy- 
sis for m x > 50 GeV/c 2 (solid red), m x = 10 GeV/c 2 (dotted 
green), m x — 7 GeV/c 2 (dash-dotted black). The optimum 
interval analysis additionally uses a S2/S1 ER discrimination 
cut. Its NR acceptance is also shown (dashed blue). 



The energy window for the WIMP search is chosen be- 
tween 4 — 30 PE, corresponding to 8.4 — 44.6 keV nr based 
on the £ e ff parametrization shown in Fig. 1. The lower 
bound was assumed to give a sufficient discrimination 
between genuine SI signals and electronic noise, whereas 
the upper bound was chosen to include most of the ex- 
pected WIMP signal. To discriminate ERs from NRs, 
the parameter log 10 (S2b/Sl) is used, with S2b being the 
sum of the bottom PMT signals. This choice is motivated 
by the more uniform signal distribution on these, giving 
a smaller uncertainty from the (x, y) S2-correction. To 
maximize the sensitivity given the homogeneously dis- 
tributed 85 Kr background, the fiducial volume was opti- 
mized on ER background data and set to 48 kg, and the 
ER rejection level was set to 99.75%. The acceptance to 
NRs below this line is calculated based on single-scatter 
NRs from 241 AmBe data and is also shown in Fig. 2. The 
Profile Likelihood approach tests the full S2/S1 space and 
does not employ this cut. 

The expected background in the WIMP search region 
is the sum of Gaussian leakage from the ER background, 
of non-Gaussian leakage, and of NRs from neutron in- 
teractions. The latter, estimated by Monte Carlo sim- 
ulations, takes into account neutron spectra and total 
production rates from (a, n) and spontaneous fission re- 
actions in the detector and shield materials, and is based 
on the measured radioactivity concentrations [8]. The 
impact of muon-induced neutrons is obtained from sim- 
ulations and contributes 70% to the total. Taking into 
account the measured trigger efficiency and the energy 



4 



threshold in the active LXe veto, the overall prediction 
is (0.31^0 n) single scatter NRs in the 100.9 days data 
sample before a S2/Sl-cut, in the energy region of in- 
terest and 48 kg fiducial mass, of which (0.11+q'q 8 ) are 
expected in the benchmark WIMP search region. 
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FIG. 3: Observed event distribution using the discrimina- 
tion parameter log 10 (S2b/Sl), flattened by subtracting the 
ER band mean, as a function of NR equivalent energy (keV nr ). 
All quality cuts, including those defined after unblinding, are 
used. Gray points indicate the NR distribution as measured 
with an 241 AmBe neutron source. The WIMP search region is 
defined by the energy window 8.4 — 44.6 keV nr (4 — 30 PE) and 
the lower bound of the software threshold S2 > 300 PE (blue 
dashed) . The optimum interval analysis additionally uses the 
99.75% rejection line from above and the 3<r contour of the 
NR distribution from below (green dotted). Three events fall 
into this WIMP search region (red circles), with (1.8 ± 0.6) 
events expected from background. 

The normalized ER band, obtained by subtracting its 
mean as inferred from 60 Co calibration data, is well 
described by a Gaussian distribution in log 10 (S2b/Sl) 
space. Gaussian leakage, dominated by the 85 Kr back- 
ground, is predicted from the number of background 
events outside the blinded WIMP search region, taking 
into account the blinding cut efficiency and the ER re- 
jection level. It is (1.14 ± 0.48) events in the benchmark 
WIMP search region, where the error is dominated by the 
statistical uncertainty in the definition of the discrimina- 
tion line. Non-Gaussian (anomalous) leakage can be due 
to double-scatter gamma events with one interaction in a 
charge insensitive region, e.g. below the cathode, and one 
in the active target volume. Such events have a lower ef- 
fective S2/S1 ratio, since only one interaction contributes 
to the S2, but both to the SI. Their contribution has 
been estimated using 60 Co calibration data, taking into 
account the different exposure compared to background 
data, and accounting for the fact that the background of 
this data set is dominated by 85 Kr which /3-decays and 
does not contribute to such event topologies. The spatial 
distribution of leakage events for background and calibra- 
tion data is similar within 10%. This is verified by Monte 



Carlo simulations and by data, selecting potential leak- 
age candidates by their SI PMT hit pattern. Anomalous 
leakage is estimated to give (0.56io"27) events, where the 
uncertainty takes into account the difference in the back- 
ground and calibration distributions, and that the leak- 
age might be overestimated because of the uncertainty 
in the 85 Kr concentration. In summary, the total back- 
ground prediction in the WIMP search region for 99.75% 
ER rejection, 100.9 days of exposure and 48 kg fiducial 
mass is (1.8 ± 0.6) events. This expectation was verified 
by unblinding the high energy sideband from 30 — 130 PE 
before unblinding the WIMP search region. The Pro- 
file Likelihood analysis employs the same data and back- 
ground assumptions to obtain the prediction for Gaus- 
sian, non-Gaussian and neutron background for every 
point in the log 10 (S2b/Sl) parameter space. 

After unblinding the pre-defined WIMP search region, 
a population of events was observed that passed the SI 
coincidence requirement only because of correlated elec- 
tronic noise that is picked up from an external 100 kHz 
source, as verified by inspection of the digitized PMT sig- 
nals. These events are mostly found below the SI analysis 
threshold, with 3 events from this population leaking into 
the WIMP search region close to the 4 PE lower bound. 
This population can be identified and rejected with a cut 
on the SI PMT coincidence level, that takes into account 
correlated pick-up noise, and by cutting on the width 
of the SI candidate. These post-unblinding cuts have a 
combined acceptance of 99.75% for NRs while removing 
the entire population of noise events. 



Radius [cm] 
10 12 




100 150 
Radius 2 [cm 2 ] 



FIG. 4: Distribution of all events (gray dots) and events be- 
low the 99.75% rejection line (black dots) in the TPC observed 
in the 8.4 — 44.6keV nr energy range during 100.9 live days. 
All cuts are used here, including the ones introduced post- 
unblinding to remove a population due to electronic noise. 
The 48 kg fiducial volume (dashed, blue) and the TPC di- 
mensions (gray) are also indicated. 

With these additional cuts, 3 events pass all quality cri- 
teria for single-scatter NRs and fall in the WIMP search 
region, see Fig. 3. This observation remains unchanged 
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for moderate variations in the definition of any of the data 
quality cuts. These events were observed on January 23, 
February 12, and June 3, at 30.2 keV nr , 34.6 keV nr , and 
12.1 keV nr , respectively. The event distribution in the 
TPC is shown in Fig. 4. Given the background expecta- 
tion of (1.8±0.6) events, the observation of 3 events does 
not constitute evidence for dark matter, as the chance 
probability of the corresponding Poisson process to re- 
sult in 3 or more events is 28%. 
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FIG. 5: Spin-independent elastic WIMP-nucleon cross-section 
a as function of WIMP mass m x . The new XENON100 limit 
at 90% CL, as derived with the Profile Likelihood method tak- 
ing into account all relevant systematic uncertainties, is shown 
as the thick (blue) line together with the expected sensitivity 
of this run (yellow/green band). The limits from XENON100 

(2010) [7], EDELWEISS (2011) [6], CDMS (2009) [5] (re- 
calculated with v CBC = 544 km/s, v = 220km/s), CDMS 

(2011) [19] and XENON10 (2011) [20] are also shown. Ex- 
pectations from CMSSM are indicated at 68% and 95% CL 
(shaded gray [21], gray contour [22]), as well as the 90% CL ar- 
eas favored by CoGeNT [23] and DAM A (no channeling) [24]. 

The statistical analysis using the Profile Likelihood 
method [17] does not yield a significant signal excess ei- 
ther, the p-value of the background-only hypothesis is 
31%. A limit on the spin-independent WIMP-nucleon 
elastic scattering cross-section a is calculated where 
WIMPs are assumed to be distributed in an isothermal 
halo with vo = 220 km/s, Galactic escape velocity v csc — 
(544jjg) km/s, and a density of p x = 0.3GeV/cm 3 . The 
SI energy resolution, governed by Poisson fluctuations of 
the PE generation in the PMTs, is taken into account. 
Uncertainties in the energy scale as indicated in Fig. 1, 
in the background expectation and in v csc are profiled 
out and incorporated into the limit. The resulting 90% 
confidence level (CL) limit is shown in Fig. 5 and has 
a minimum a = 7.0 x 10~ 45 cm 2 at a WIMP mass of 
m x = 50 GeV/c 2 . The impact of £ e ff data below 3keV nr 
is negligible at m x = 10GeV/c 2 . The sensitivity is the 
expected limit in absence of a signal above background 
and is also shown in Fig. 5. Due to the presence of 
two events around 30keV nr , the limit at higher m x is 



weaker than expected. Within the systematic differences 
of the methods, this limit is consistent with the one from 
the optimum interval analysis, which calculates the limit 
based only on events in the WIMP search region. Its 
acceptance-corrected exposure, weighted with the spec- 
trum of a m x = 100GeV/c 2 WIMP, is 1471kg x days. 
This result excludes a large fraction of previously unex- 
plored WIMP parameter space, and cuts into the region 
where supersymmetric WIMP dark matter is accessible 
by the LHC [21]. Moreover, the new result challenges 
the interpretation of the DAMA [24] and CoGeNT [23] 
results as being due to light mass WIMPs. 
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